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Abstract 

An improved formulation of the one-step model of photoemission from crystal surfaces 
is proposed which overcomes different limitations of the original theory. Considering the 
results of an electronic-structure calculation, the electronic (one-particle) potential and the 
(many-body) self-energy, as given quantities, we derive explicit expressions for the dipole 
transition-matrix elements. The theory is formulated within a spin-polarized, relativistic 
framework for general nonspherical and space-filling one-particle potentials and general 
nonlocal, complex and energy- dependent self-energies. It applies to semi-infinite lattices 
with perfect lateral translational invariance and arbitrary number of atoms per unit cell. 

PACS: 79.60.-i, 73.90.+f 
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1 Introduction 



The spectrum of one-particle excitations of a metallic system of correlated electrons in a 
solid is a fundamental question in condensed-matter physics. Experimentally, the inter- 
esting valence-band region around the Fermi energy is accessible by means of ultraviolet 
photoemission spectroscopy (PES) Q and inverse photoemission spectroscopy (IPE) 0. 

The theoretical understanding of the excitation spectrum poses a long-standing and 
not yet generally solved problem. Within the independent-electron approximation the 
spectrum is simply given in terms of the one-particle (Hartree-Fock) eigenenergies of the 
Hamiltonian. Analogously, it is widely accepted to interpret a measured photoemission 
spectrum by referring to the results of band-structure calculations that are based on den- 
sity functional theory (DFT) and the local density approximation (LDA) [|3], ||. Despite 
a sometimes convincing success in practice @, 0], such an interpretation is questionable 
since there is actually no known correspondence between the Kohn-Sham eigenenergies 
and the one-particle excitations of the system ||, 0]. For an in principle correct descrip- 
tion of the excitation energies, the local LDA exchange-correlation potential has to be 
supplemented by the nonlocal, complex and energy-dependent self-energy which leads to 
the Dyson equation ||, || instead of the Schrodinger-type equation in the Kohn-Sham 
scheme. 

Provided that the self-energy is known, one can deduce a PES/IPE "raw spectrum" 
from the solution of the Dyson equation. To achieve a reliable interpretation of exper- 
iments, however, it is inevitable to deal with so-called "secondary effects" which con- 
siderably modify and distort the raw spectrum. Above all, the wave-vector and energy 
dependence of the transition-matrix elements has to be accounted for. These dependencies 
are known to be important and actually cannot be neglected. They result from strong 
multiple-scattering processes which dominate the electron dynamics in the low-energy 



regime of typically 1-100 eV ||10|| . The transition-matrix elements also include the effects 
of selection rules which are not accounted for in the raw spectrum. Strictly speaking, 
it can be stated that the main task of a theory of photoemission is to close the gap be- 
tween the raw spectrum obtained by (many-body) electronic-structure calculations and 
the experiment. 

It turns out, however, that the calculation of the transition-matrix elements is by 
no means a trivial task. Probably, the most successful theoretical approach is the so- 
called one-step model of photoemission as originally proposed by Pendry and co-workers 
]T(]| , |TT| , |r2|] . A review on the recent developments and refinements |T3| of the approach 
can be found in Ref. 0. The main idea of the one-step model is to describe the actual 
excitation process, the transport of the photoelectron to the crystal surface as well as the 
escape into the vacuum |14| as a single quantum-mechanically coherent process including 
all multiple-scattering events. 

The main disadvantage of the conventional formulation of the one-step model consists 
in the fact that it is intrinsically based on a local potential. This is sufficient if the 
calculation starts from the self-consistent LDA potential. The inclusion of the nonlocal 
self-energy which is needed for an in principle correct description of the one-particle 



2 



excitations, however, is not possible. In a preceding study |15j we were able to show that 
one can overcome this difficulty and that the nonlocal self-energy term can be included 
within an alternative formulation of the one-step model. 

The purpose of the present paper is to generalize the results of Ref. [ITS] in different 
respects: (i) The two-component formalism is replaced by a four-component, relativistic 
framework necessary for the study of high-Z materials or for spin-polarized photoemission 
from nonmagnetic samples excited by circularly polarized radiation, for example, (ii) For 
ferromagnetic systems, the exchange splitting and the spin-orbit splitting are treated on 
equal footing. This allows to investigate e. g. magnetic dichroic effects, (iii) The theory 
is no longer based on the muffin-tin approximation for the input LDA potential. Instead, 
the general case of space-filling, nonspherical potentials is considered which becomes im- 
portant for more open crystal structures, (iv) To cover the case of complex geometries 
(ordered compounds, multilayers), we also generalize to more than one atom per unit cell, 
(v) The Korringa-Kohn-Rostocker (KKR) multiple-scattering formalism [16, 17 is used 
for both, the final as well as the initial states. This internal consistency of the formalism 



represents another improvement compared with Ref. |T5| where a muffin-tin-orbitals basis 
was used with respect to the initial states. 

As in Ref. |15[ we keep the basic structure of the one-step model. Starting from 



Pendry's formula for the photocurrent ||11|| , analytical expressions for the matrix elements 
are derived referring to a crystal surface with perfect lateral translational invariance. 
For the solution of the (atomic) Dirac equation we use the phase-functional ansatz of 



Calogero generalized to the nonspherical case [JT8J, [L9|j . The initial state is treated within 
the relativistic version of KKR theory for space-filling potentials |2(J [H| but adapted 
for a slab geometry, and the final state is constructed using a full-potential relativistic 
layer-KKR method [^, ^3|. The new formalism is developed up to the point where the 
numerical evaluation has to start. 

The LDA potential and the self-energy are assumed to be given (input) quantities 
that must be obtained from a preceding electronic-structure calculation. The related 
many-body problem is beyond the scope of the present paper. In particular, one must 
be aware of a possible double counting of interactions, once on a mean-field level in the 
LDA and once explicitly in the self-energy. While there are pragmatic ways to circumvent 
this problem p^ , f25f , it generally remains unsolved. Whether or not the self-energy 
is calculated within a relativistic framework is not decisive for the applicability of the 
presented approach. A two-component, nonrelativistic self-energy as it is usually provided 
(see e. g. p4|, [26|, ^7j) can be used as an input quantity in the relativistic formulation. 



2 One-step model of photoemission 

Photoemission (PES) and inverse photoemission (IPE) are complemental spectroscopies. 
We concentrate on PES in the following since IPE can simply be treated analogously 
by taking into account geometrical factors that regard the respective experimental setups 
p8| . We start our considerations by a discussion of Pendry's formula for the photocurrent 
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which defines the one-step model of PES ||11||: 



J PES oc Im (e^kiilG+AG+A^le^kii) . (1) 

The expression can be derived from Fermi's golden rule for the transition probability 
per unit time §|. Consequently, J PES denotes the elastic part of the photocurrent. Vertex 
renormalizations are neglected. This excludes inelastic energy losses and corresponding 



quantum-mechanical interference terms || |IT| |29|. Furthermore, the interaction of the 
outgoing photoelectron with the rest system is not taken into account. This "sudden 
approximation" is expected to be justified for not too small photon energies. 

We consider an energy-, angle- and spin- resolved photoemission experiment. The state 
of the photoelectron at the detector is written as |e/, kn), where kii is the component of 
the wave vector parallel to the surface, and e/ is the kinetic energy of the photoelec- 
tron. The spin state of the photoelectron is implicit in |e/, ky) which is understood as 
a four-component Dirac spinor. The advanced Green function G 2 in Eq. ([!]) charac- 
terizes the scattering properties of the material at the final-state energy E 2 = e/. Via 
\f) = G2 \tf, k||) all multiple-scattering corrections are formally included. Using standard 
Korringa-Kohn-Rostocker (KKR) multiple scattering techniques p2| , we can calculate the 



final state |/) as a (time- reversed) relativistic LEED state (see Sec. (|). 

As far as concerns the final state, many-body effects are included only phenomenolog- 
ically in the LEED calculation, i. e. by using a parametrized, weakly energy- dependent 
and complex inner potential V (E 2 ) = V 0r (E 2 ) + iV i(E 2 ) as usual [[TTJ. This generalized 
inner potential also includes the (imaginary) optical potential, which takes into account 
inelastic corrections to the elastic photocurrent |§ as well as the actual (real) inner poten- 
tial, which serves as a reference energy inside the solid with respect to the vacuum level 
fjQfl . Due to the finite imaginary part iVai(E 2 ), the flux of elastically scattered electrons 
is permanently reduced, and thus the amplitude of the high-energy wave field |/) can be 
neglected beyond a finite distance from the surface. It is thus sufficient to restrict oneself 
to a slab geometry in a practical computation. 

A in Eq. ([!]) is the dipole operator in the electric dipole approximation which is well 
justified in the visible and ultraviolet spectral range. It mediates the coupling of the 
high-energy final state with the low-energy initial states. For the relativistic, possibly 
ferromagnetic case and for general space-filling potentials, a convenient form of the dipole 
operator is given in Sec. ||. 

The "low-energy" propagator Gf in Eq. ([]]), i. e. the one-electron retarded Green 
function for the initial state in the operator representation, yields the "raw spectrum" . It 
is directly related to the "bare" photocurrent and thereby represents the central physical 
quantity within the one-step model. Gf = Gf(Ei) is to be evaluated at the initial-state 
energy Ei = e/ — uo — fiQ, where u is the photon energy (/i stands for the chemical 
potential). 

In the framework of the conventional one-step model of Pendry and co-workers [|H], 
0] the initial-state Green function Gf is determined for f, J, electrons moving in 
an (effective) one-particle potential VldaM provided by DFT-LDA. As usual, VldaM 
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consists of the external core potential, the Hartree contribution as well as the exchange- 
correlation potential. In the relativistic generalization of DFT [[H], one has to consider 
the (one-particle) Dirac Hamiltonian (h = m = e = l,c = 137.036): 

h LDA (r) = -tcctV + (3c 2 - c 2 + V^a^) + (3(rB LDA (r) , (2) 

where Vlda(i") denotes the (effective) spin-independent potential, and Blda(i") is the 
(effective) magnetic field. They are given as j33| : 



Wr) = -(^iIdaW + V^ DA (r)) (3) 

and 

B LDA (r) = ^ L T DA (r) - Vjf DA (r)) b . (4) 

The constant unit vector b determines the spatial direction of the (uniform) magnetization 
as well as the spin quantization axis. f3 denotes the usual 4x4 Dirac matrix with the 
nonzero diagonal elements /3n = P22 = 1 and ^33 = /J44 = —1, and the vector a is given 
by its components otk = o x ® (k = x, y, z) in terms of the 2x2 Pauli-matrices o^. 

Within the DFT ground-state formalism, VldaM as well as B LDA (r) are local func- 
tions. On the other hand, it is well known that for an in principle exact description of the 
one-particle excitations one has to consider the Dyson equation for the Green function 
[|. [|. This includes the nonlocal, complex and energy-dependent (retarded) self-energy 
ETi(r, r ' As the LDA potential, the self-energy must be assumed to be a given quan- 
tity for the photoemission theory. Therewith, we can construct a generalized potential, 

U(r,r',E) = 6{r - r') (V LDA (r) + /3<rB LDA (r)) 

+ V{r,r',E) + /3aB(r,r',E) , (5) 

where the nonlocal contributions V and B are defined as: 

V(r,r',E) = ~(E T (r,r',£) + E*(r, r', E)) (6) 

and 

B(r,r',£) = ~(E^r, r', E) - E^r, r', E)) b . (7) 

The initial state is described by a 4 x 4 Green matrix Gi(r, r', E\). It can be obtained 
as the solution of the Dyson equation which can be written as: 

[E l + f i -h hBA (r)]Gt(r,r',E 1 ) J (V(r, r", E x ) + /3<rB(r, r", E 1 )) 

* G+(r", r', Ejdr" = 8(r - r')l . (8) 

A direct solution of the Dyson equation in real-space representation turns out to be 



inconvenient. Analogously to Ref. [15], we therefore turn to a matrix representation and 
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choose the eigenspinors of the LDA Hamiltonian (0) as basis states. Assuming perfect lat- 
eral translational symmetry, the parallel component of the wave vector is a good quantum 
number, and the LDA eigenvalue problem reads: 

/iLDAKqy) = e n (q||)|n, qy) . (9) 

Here qy is a vector of the first two-dimensional Brillouin zone, and e„(qy) is the two- 
dimensional band structure. Using the eigenspinor basis, the Dyson equation for the 
initial-state Green function can be written: 

]T ({Ex + mo - en(q||))<W - V nn .{Ex, qy) - B nn ,{E x , qy)) G%l„{E 1: qy) = 5 nn „ . (10) 



Here we have introduced the matrix representation of the nonlocal terms V and B: 

Kn'(#i,q||) = (n,q|||V(£i)K,qy), 

B nn ,(Ei,q||) = (n.qnl^B^OI^qii). (11) 

Lateral translational symmetry requires V and B to be diagonal with respect to qy . 
In the eigenspinor basis Pendry's formula reads: 

Joclm ^M n ( e/ ,k||) G^(^i.qil) K'( e /> k ll) > (12) 

nn' 

where 

M ri ,(e / ,k|,) = (e / ,k|||G' 2 + A|n,q||) (13) 
is the matrix element of the dipole operator between the final state |/) = G2 |e/, ky) and 



the LDA eigenspinor \n, qn). In Eq. (12) qy is fixed by translational symmetry: A nonzero 



matrix element fll3|) requires qn = ky + gy, where gii is a reciprocal lattice vector (15 . 
Furthermore, ky is given by the photoelectron energy at the detector and by the emission 
angles. This implies that gy and thus qy = qy(ky) are uniquely determined by demanding 
qy to lie within the first two-dimensional Brillouin zone. 

To work out the differences with respect to the original formulation, let us briefly 
discuss the case S = 0. Energy and momentum conservation then implies that the final 
state \ f) can couple via A to a single "initial" spinor |n, qy) only: Its eigenenergy is given 
by e n (qy) — €f — u, and qy = ky + gy. This implies that we can set 

G t = K p 7 — 7 — r x ~ - n+ (n, qy I (14) 

Ex - (e n (qy) - fi ) + i0 + 

in Eq. ([[]). I. e. the conventional one-step model is characterized by a one-pole structure of 
the low-energy (hole) propagator. It is because of this (implicit) one-pole structure that 
the original formulation does not allow to include a general nonlocal self-energy term. 
(Essentially the same argument applies if the infinitesimal i0 + is replaced by a small 
finite imaginary constant as is oftenly done.) 
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On the contrary, for the general case H = V + By^0we have: 



G+ = $>,q||) 



11 

E 1 -(e(q||)- / i )-S(^i,q||) 



(n',qi||. (15) 



The imaginary part of the self-energy causes an energy, wave- vector and band dependent 
broadening of the initial state. The hole acquires a finite lifetime (except for E\ = 0, i. e. 
for €f — uj = no where the self-energy is Hermitian). The energy dependence of E may 
introduce a band-narrowing effect, an enhancement of the effective electron mass at the 
Fermi edge and, for a strongly correlated system, may give rise to satellite features in the 
spectrum. All this implies that for a given photon energy uj the final state |/) necessarily 
couples to different initial states |n, qy). Furthermore, since £ is generally nondiagonal in 
the band index n, interference terms ?i/n' have to be considered. 

The number N n of possible values for the band index n in Eqs. (|10D and ([12]) is given 
by N n = Nj_N^N^ where Nj_ is the number of layers in the slab and is the number of 
atoms in the two-dimensional unit cell. is determined by the maximum order in the 
spin-angular momentum expansion (see next section) that is necessary for convergence. 
N± can be assumed to be finite since the damping of the final-state wave field implies 
that contributions to the photocurrent are negligibly small beyond a finite distance from 
the surface. In fact, a finite Nj_ is decisive for a numerical solution of the Dyson equation 
(|10|) by matrix inversion. One inversion is necessary to get one (e/,kii) point in the PES 
spectrum. 



Formulas (0), ([12] ) and (|13|) generalize our approach of Ref. ||15| . The latter is obtained 
in the nonrelativistic approximation. On the other hand, in the limit S^(r, r', E) = the 
initial state is treated as in the original (relativistic) formulation of the one-step model 



The actual and remaining task, however, consists in the according recalculation of 
the transition- matrix elements (0). This is done in the following. We thus consider the 
initial states \n, qy) (Sec. ||), the final state G^|e/,k||) (Sec. f|), and the dipole operator 
A (Sec. |5|) to derive the final expression for the matrix elements in Sec. []. 



3 The initial states 

For the evaluation of the matrix elements, we need a one-center expansion of the LDA 
eigenstates \n, qii). This will finally allow for a separation of the matrix elements into 
radial and angular parts. Such a one-center expansion is available within the relativistic 



version of the KKR approach for arbitrary space-filling potentials |17 . 

We consider the slab to be built up from layers i± = 1, N± parallel to the surface. 
The two-dimensional unit cells within a layer i± are labeled by an index i\\ and the atoms 
within the unit cells by an index The position vector of a particular atom in the 
semi-infinite lattice % = (i\\,i±,iA) is then given by: Rj = Rj.. + R i± + R iA , where Rj,. is 
a vector of the two-dimensional lattice, Rj x denotes the local origin of layer and R« A 
the position vector of the «A-th atom with respect to the local origin of the unit cell. The 
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material space is decomposed into (three-dimensional) polyhedra (atomic cells) Qi with 
one atom i at the center of each such that the LDA potential Vlda(f) (and analogously the 
effective magnetic field Blda(f)) can be written as a sum over cell potentials Vi(r — Rj) 
that vanish outside Qf 

V LDA (r)=J2V l (r-R 2 ). (16) 

i 

A cell Qi is circumscribed by a bounding sphere Si of radius Ri. Due to translational 
symmetry the potential within the cell Qi or within the sphere Si only depends on i± and 
za- We have: V*(r) = V^ x i A (r). The same notation also applies for the effective magnetic 
field etc. 

For the practical calculation, slabs consisting of N± layers each are arranged in a 
super-cell geometry with a sufficiently large distance in between. This formally restores 
full three-dimensional periodicity and thus allows the application of the conventional KKR 
method. It can be shown J2l|] that the wave function \&W(r) = (r|n, qy) of the initial 
spinor |n, qy) can be expanded in locally exact basis spinors Q iKfl (E,r) = $ ixiAKfl (E,r): 

*W(r) = J2 A^(e n (q||), q„) ^(e n (q ]{ ), r - R) . (17) 

Kfl 

The expansion converges within each cell fli and for each e n (q||) [21p. The basis spinors 
satisfy 

(E + icotV - (3c 2 + c 2 - K(r) - PtrBiir)) & iKll (E, r) = (18) 
for all r. They differ in their behavior close to the center of For |r| i— > 0: 

^(£,r) = J£(*,r), (19) 

where 

J " {K r) = ( ix^{M)rii kr ) ) ' (20) 



with the spherical Bessel function ji(kr) and k, S K , and I defined by k = y2E + E 2 /c 2 } 
S K = «/|«|, and I — I — S K . k and /i are the relativistic spin- angular momentum indices 



according to Rose ]34}1 . The spin-angular functions = C^Yf 1 ~ s (r )x s are given 

in the usual way [34 in terms of Clebsch-Gordan coefficients C KfJ-s , spherical harmonics 
Y t m and Pauli spinors X1/2 = (1,0)*, X-1/2 = (0, l) f . 

Eq. (p~8f) is a well-defined (effectively atomic) problem which together with fll9l) uniquely 
determines the &i Kll (E,r). It can be solved using the (nonspherical) phase-functional 
ansatz of Calogero generalized to the relativistic case |18|, [19T| . We separate radial and 
angular parts and write: 

= Y,(4{K*)Ci.^>n,,{E,r)-N£{kir)S iK ^,^{E,r)) , (21) 
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where N%(k, r) is defined analogous to Eq. ( p0|) replacing ji(kr) by the spherical Neumann 
functions ni(kr). The superscripts u and / refer to the upper and lower components of 
the four- component spinor, respectively. Inserting into fll8|) eventually yields the coupled 
channel equations jnj] for the coefficient matrices Cj and 5j. The coupled channel equa- 
tions may be solved by outward integration from the origin to the radius of the bounding 
sphere _R, for any effective potential that is less singular than r~ 2 (p = k(E + 2c 2 ) /c): 

C^ v (£,r) = 5 KK ,5^-p f r 2 dr [ dr N£{k, r)t (^(r) + /fcrB^r)) $^,(£7, r) , 

JO J(4-n-) 

S^ KV (£7,r) = -p I' r 2 dr [ dr J%(k, r)t (^(r) + /^(r)) $« v (£, r) . (22) 

JO J(47r) 

Expanding V^(r) and Bj(r) in spherical harmonics, the angular integrations can be done 
analytically. Outside the bounding sphere where V^(r) = Bj(r) = 0, the coefficient 
matrices C iK ^^(E,r) and S in ^ K >^{E,r) are constant. 

Using the result (p2|) in Eq. and inserting into Eq. (|17D yields the desired one- 
center expansion of the initial states. The coefficients A iKAt (e n (q||), qn) in ( fPTD can be 
obtained from the KKR secular equation which for a slab geometry reads: 

— X! ^i±«A«WxV>«"M"(-^> ^11 ) ^Vk'V'kV^' I Aix'iA'K'n'(E, <l\\) = . (23) 

£> are the usual (three-dimensional) KKR structure constants (cf. Ref. |]21| ) for the super- 
cell geometry (each super-cell consists of iVj_ x TYa atoms). 



4 The final state 

The final state (r) = (r|(?2 |e/, ky) is an eigenspinor of /ilda with eigenenergy _E 2 = e/ 
and could thus be constructed in the same way as the initial states. However, for an 
appropriate description of the photoemission process we must ensure the correct asymp- 
totic behavior of ^ (r) beyond the crystal surface, i. e. a single outgoing plane wave 
characterized by e/ and ky . Furthermore, the damping of the final state due to the imag- 
inary part of the inner potential iVo\(E2) must be taken into account. We thus construct 
the final state within SPLEED theory considering a single plane wave |e/, kn) advanc- 
ing onto the crystal surface. Using the standard layer-KKR method [TH] generalized for 



the relativistic, full-potential case (cf. e. g. Ref. 0), we first obtain the SPLEED state 
U^if (r). The final state is then given as the time- reversed SPLEED state (U = —ia y K 
is the relativistic time inversion). 

The scattering properties of an atomic cell fli can be determined from the solution of 



the coupled channel equations fl22|) for the final-state energy e/. From C iKflK ' ^(E = €f,r 
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Ri) and S iK ^ K '^(E — €f,r — R{) the atomic scattering matrix T iKllti i^ is constructed as: 

where the new coefficient matrices U{ and are given by: 

Ui K ^ K '^i = Ci^^i^i (e^, R{) + i Si KfJLK i^i (e^, , (25) 
Kk/ukV = Ci^'n^ef, Ri) — i Si KllK i^{ef,Ri) . (26) 

The matrix Tj = rj ± j A together with the crystal geometry determine the scattering matrix 
M i± for a single layer i± |T|: 

"^gz K/i K »^» i A i A > 

* r, xiAl ^ v «(i - x)ih^^,^/c^ s Y{?- s '(^y^^ . (27) 

k^ denotes the wave vector and r = ± defines the direction of a plane wave (incoming or 
outgoing) with respect to the layer i±. The parallel component of kg differs from ky by 
a two-dimensional reciprocal lattice vector g while the perpendicular component is fixed 
by the plane wave energy ej. Multiple scattering processes within the layer i± are taken 
into account via the X matrix of layer-KKR theory |TD|, [TB[ : 



* G ? i/Lt-s,/"/i"-s(R-i A ' - R-ia — Rin)CK'v"sr i _ 1 _ iA ' )K; // At // /t / M / . (28) 

The prime on the summation over i\\ indicates that the term where Rj,, = and at the same 
time i/\ = i/\ is omitted. An explicit expression for the lattice sum Gi m ^i m i (Rj A '— Rj A — Rj.. ) 
can be found in Ref. [TQ . 

Therewith, the scattering properties of all layers are known. We consider a plane wave 
|e/,k||) advancing onto the crystal from the vacuum side. From the scattering matrices 
(p7|) the coefficients w[ ±gs of an expansion into plane waves |g, r, s) in front of all layers 
i±_ can easily be found using standard recursive layer-by-layer schemes [H| 2^]. Time 
reversal then yields the final-state wave function \&j^(r). 

Within each atomic cell Qj the final state may be expanded, 



< } (r) = E k|])* *S( £ /. r " R *) > ( 29 ) 



in the locally exact basis of time-reversed ( "T" ) phase functions: 



o^)=e( W!'t ) • (30) 

1 (r) 0L VreM (^,r)* 
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Here X« (t) * s a time-reversed spin-angular function: 

X^ iT) (r) =UxZ{*) = J2(-2s)C K ^ s Yf +s (r)* Xs . (31) 

s 

The expansion coefficients in Q25| ) are derived within SPLEED theory as usual. Starting 
from the bare coefficients, 



<w = EE tt, '(-^)(-)^ s C fi y s ^ A 



«V g 



/../ ^ 



[< gs ^'(Q) e^ + ' R - + u- ±ss Yr»\k s ) e°*-"*] , (32) 
and correcting for intra-layer multiple scattering, we obtain: 

Ai ± i AKfJ , — ^ ^2 ^2 ^2 ^jlLkMkV (1 _ ^)^x 1 ,^A«'M^W«''M''^ i A^«''M' / «'''M'''A■ x ^ A '«'V'' ' 



The one-center expansion (|29|) is formally very similar to the one-center expansion 
(O) of the initial states. The main difference, however, is the calculation of the expansion 
coefficients. For the final state they are obtained from (full-potential) SPLEED theory 
instead of solving the KKR secular equation. Thereby, we account for the final-state 
damping effects and the correct asymptotics corresponding to the experimental situation. 



5 The dipole operator 

In the relativistic theory the dipole interaction of an electron with the electromagnetic 
field is given by the dipole operator A = — aA where A is the spatially constant vector 
potential inside the crystal. In a matrix element (^|A|^) between eigenspinors 
and \^>i) of the Dirac Hamiltonian (0) with energies Ef and Ei, respectively, A can be 
written as: 

A(r) = E fi (Ac V + ^qA ) V LDA (r) + E fi (A V) /3<xB LDA (r) + E fi ^(3A x <xB LDA (r) , 

(34) 

with Efi = —2ic/[(Ef + c 2 ) 2 — (Ei + c 2 ) 2 ]. The expression is derived by making use of 
commutator and anticommutator rules analogously to the nonrelativistic case in Ref . |15|] . 
Using the decomposition of the LDA potential and the effective magnetic field (|T6|) , we 
immediately get: 

A(r)=]TA ? ;(r), (35) 

i 

where Aj(r) vanishes outside flj. In the following we calculate the dipole operator for a 
given atomic cell. Distinguishing between contributions due to V{(r) and due to B.;(r), 
the 4x4 matrix in (Bl) can be written more explicitly as: 



A^(D A;'!/(r) 
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with the following 2x2 matrices: 



A( uu ) /„ 
v,i ( r 



A W;(r) ]l (2x2) 
A£ ) W = Ag ) (r) 



Ag?(r) = (A V)( < rB i (r)) - ^A x B,(r) 



A2J(r) 







Ag } (r) 
AiS(r) 



Ag?(r) 
Ag(r) 



0-A o Vj(r) 



(37) 







-A 



H, 



As for the initial states and the final state we need a separation into radial and angular 
parts. For this purpose the cell potential Vi(r) inside the bounding sphere is expanded 
into spherical harmonics, Vi(r) = Ezm^mM 5 ?™^) ( and B *( r ) = b T,i m B i i m {r)Y l m (r)). 
The expansion has to be done in a way that guarantees VldaC 1 ") = B L DA( r ) = in the 
segments between the bounding sphere and the atomic cell (Si — 0$). For the two terms 
in (|37|) involving the gradient the expansion needs to be considered in detail: 



AoW 4 (r) = EE^W^( f ) > 

a=l Im 

(A V)(<xB 4 (r)) = EE5£(r)(<T-b)^(f) 

a=l lm 



A straightforward calculation yields 



ilm 



BZ(r) 



for the radial parts vj^l and due to Vi(r) and Bj(r), respectively, and 



' d_ i + r 



Bum(r) ■ 



-K---] V Um (r) 
or r 



B Sl(r) = [l - l A B ilm (r) 



W^(r) 



m'=l 



E n m '(^,) X £L Y^t m '(i) 



m'=-l 
m'=l 



m'=-l 



with 



Y (l) 



(2) 




(/ — m)(l + m) 
\l (2Z-1)(2Z + 1)\ 



2/ - m'(m + m'(l + 1)) 
2(Z + m'm) 




(I + 1 - w)(Z + 1 + m) 



4tt\ (2/ + l)(2/ + 3) \ 2(/ + l -m'm 



2(Z + 1) + m'(m - m'l) 



(38) 
(39) 

(40) 
(41) 



(42) 



(43) 



for the angular parts wj^ (a = 1, 2). is the unit vector in the direction of the vector 
potential A . 
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6 The transition-matrix elements 



Now we are in the position to calculate the transition-matrix elements 

M n ( e/ ,k||) = ( e/ ,k|||G+A|n,q||) = J d 3 r (r)tA(r)*<J>(r) . (44) 



Inserting the one-center expansions for the initial and the final states ( |17D and (p9| 
respectively, and using the decomposition (|35|) of the dipole operator yields: 



M n (e/,k|| 



fan)' mi 

N \\H E A ixiAKfM (e f ,k l} ) M,_ 



/(e / ,e n (q||)) A i± 

iak' a*' ( e n(q||) j qn 



where iVn is the number of two-dimensional unit cells per layer and 



(45) 



(46) 



is the atomic matrix element of the dipole operator between the phase functions at the 
respective eigenenergies. The integration extends over the volume of the polyhedron f^. 
Since V^(r) = Bj(r) = and thus Aj(r) = in the segments between fli and the bounding 
sphere Si, the integration may be performed over Si as well which is more convenient. 
In the following we suppress the atomic index i. Using (^) we get eight dipole matrix 
elements, 

M K , K y(E f , Ei) = E fi E / d'r *$p{E f , r)t (A^(r) + A { f\r)) r) , (47) 

dd> Js 

each of which can be separated into radial and angular parts: 



M K/J , K i^(Ef, Ei) = E f i E El 

k"h"k'"h"' Im I 

E 



a=l 



„//„/// 



W («J) 



Vi 



I III 



At* A* 



+ E 



(m) 



5 



(a) 



/m 



4 M 



E rcy 



(«) 



a=l 



a=l 



7 ^ 



vi a) 



E b\ 



0=1 



,(«) 



Irn 



U ->(«) 



,4' 



//„/// 



jm'V" 



<xA x bY^ 



<xA x bY, m 



(48 
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with I = (KfiK 1 'fit 'k" >" '«'" >'") and: 



Ft 



r dr 



K" IX" K/J, 



r,E f ) f(r) cf ) i,,„,, l {r,E l ) 



AM 



f 



(4tt) 



(49) 
(50) 



for d = u,l and d' = u, I. 

This completes our formalism. From the different contributions (|48| ) we can calculate 
the atomic matrix element (|46D for each layer i± and each atom i& in the two-dimensional 
unit cell. Inserting into fl4"5]), one obtains the full transition-matrix element ([13]) which 
enters the expression ([12]) for the PES intensity. For a single point (ey, kn) in the spectrum, 
it is necessary to determine the matrix elements at the energies Ej = 6f and E{ = e n (qii) 
for each n (qy is fixed) where there is a nonnegligible contribution of the imaginary part 
of the initial-state Green function G nn '(Ei, qn). 

The different types of radial matrix elements ( |4*9"D result from the different possible 
combinations of the upper and lower components of the Dirac spinor. In the nonrelativistic 
case only the d = d! = u term survives. All one-dimensional integrals are well defined 
since for r i— > there is a r~ 2 singularity of f(r) at most. This can be seen from Eq. flilD . 
Furthermore, the radial parts of the phase functions are regular at r = [Eq. (|i~9|)1. 

For the paramagnetic case b = there are two types of angular matrix elements 
(0). The first in Eq. 



wj^ , is well known from the nonrelativistic theory |23 1 . For 



in (|48|)1 only the term for 



spherically symmetric potentials V^(r) = V m (r)Y °(r) [I = m 
a = 2 yields a nonvanishing contribution. In this case the angular matrix element further 
reduces to the muffin-tin form fll5" |. The second type of matrix elements, ~ crAol^ m , 
represents a relativistic correction [p3[| . In the ferromagnetic case two new types of matrix 
elements, ~ crbW^ , and ~ crAo x bY™, occur. All angular matrix elements are simply 
related to usual Gaunt coefficients and can thus be calculated analytically. 



7 Summary 

In this paper we have described in detail a new formalism to evaluate Pendry's formula 
for the one-step model of (inverse) photoemission spectroscopy. The initial-state Green 
function represents the central physical quantity of the one-step model. Compared with 
previous work, we propose a different treatment of the initial state and determine the 
Green function from the Dyson equation rather than within the DFT-LDA ground-state 
theory. Generally, this is a necessary condition for a correct description of one-particle 
excitations. Given the electronic self-energy, our formalism allows to include the corre- 
sponding many-body effects such as temperature dependencies, quasi-particle damping, 
band narrowing, satellites, etc. in the one-step description of the photoemission process. 



Within the framework of the original theory [10, 11, O] and its subsequent improvements 



|TB| , one is unable to deal with an in general nonlocal self-energy correction since this is 



incompatible with the KKR multiple-scattering formalism on which the original theory 
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is solely based |fLq| . Our main idea is thus to disentangle the calculation of the Green 
function from the calculation of the dipole matrix elements. This is achieved by consider- 
ing the low-energy LDA eigenstates merely as basis states to set up the Dyson equation. 
Solving the Dyson equation to get the PES raw spectrum is then independent from the 
calculation of the actual dipole matrix element. Note that it is necessary to work within 
an eigenstate basis of the LDA Hamiltonian since this leads via Eq. (|34|) to the important 
spatial decomposition fl35|) of the dipole operator which eventually allows for a convenient 
calculation of the transition-matrix elements. 

The present paper has shown that the alternative evaluation of the one-step model 
can formally be performed consistently for very general cases, namely for complex geome- 
tries with more than one atom per unit cell, for general space-filling nonspherical LDA 
potentials, for situations which require a fully relativistic treatment, and for the ferro- 
magnetic case, in particular when exchange and spin-orbit splitting must be treated on 
equal footing. 

It goes without saying that the presented formalism can also be applied without taking 
into account the self-energy correction. Even in the case E = it should be advantageous 
since it is more transparent and more simple compared with the corresponding (full- 
potential, relativistic) theory |23| based on the original concept fllf : (i) Simple one- 
dimensional integrals (|49D have to be performed instead of two-dimensional integrations 
P5fl for the calculation of the radial matrix elements. Furthermore, their number is 
reduced considerably, (ii) Only the regular instead of both, the regular and the irregular 
solution of the atomic Dirac equation is needed [see Eq. (0)]. (iii) There are no additional 



intra-layer and inter- layer contributions as in the original theory fl23| . The partitioning 



of the photocurrent into atomic, intra-layer and inter-layer contributions in the original 
formulation is only due to formal reasons, and a direct physical interpretation is difficult, 
(iv) Contrary, the separate calculation of the initial-state Green function and the matrix 
elements allows to distinguish clearly between the raw spectrum and its modifications due 
to secondary effects. 

The ferromagnetic Gd(0001) surface represents a prototype system for the application 
of the theory. Since Z = 64 for Gd, a fully relativistic evaluation of the one-step model is 
meaningful even if it is based on a potential and self-energy input from scalar-relativistic 
LDA calculations |25| and from two-component many-body theory |[25|| , respectively. Re- 
cent many-body calculations show that striking temperature-dependent correlation effects 
should be observable in the (5d, 6s) conduction band as a consequence of intra-atomic ex- 
change coupling to the subsystem of localized and ferromagnetically ordered 4f moments 
P5] , |55| . The interpretation of several PES/IPE spectra from Gd (cf. e. g. Ref. [ft5]] ) is 
controversial up to now and may be resolved by supplementing the electronic-structure 
calculations with a reliable treatment of the secondary effects. Calculations based on the 
presented formulation of the one-step model are intended for the future. 
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